56 research outputs found
A low Mach enthalpy method to model non-isothermal gas-liquid-solid flows with melting and solidification
Modeling phase change problems numerically is vital for understanding many
natural (e.g., ice formation, steam generation) and engineering processes
(e.g., casting, welding, additive manufacturing). Almost all phase change
materials (PCMs) exhibit density/volume changes during melting, solidification,
boiling, or condensation, causing additional fluid flow during this transition.
Most numerical works consider only two phase flows (either solid-liquid or
liquid-gas) for modeling phase change phenomena and some also neglect
volume/density change of PCMs in the models. This paper presents a novel low
Mach enthalpy method for simulating solidification and melting problems with
variable thermophysical properties, including density. Additionally, this
formulation allows coupling a solid-liquid PCM with a gas phase in order to
simulate the free surface dynamics of PCMs undergoing melting and
solidification. We revisit the two-phase Stefan problem involving a density
jump between two material phases. We propose a possible means to include the
kinetic energy jump in the Stefan condition while still allowing for an
analytical solution. The new low Mach enthalpy method is validated against
analytical solutions for a PCM undergoing a large density change during its
phase transition. Additionally, a few simple sanity checks are proposed to
benchmark computational fluid dynamics (CFD) algorithms that aim to capture the
volume change effects of PCMs
A moving control volume approach to computing hydrodynamic forces and torques on immersed bodies
We present a moving control volume (CV) approach to computing hydrodynamic
forces and torques on complex geometries. The method requires surface and
volumetric integrals over a simple and regular Cartesian box that moves with an
arbitrary velocity to enclose the body at all times. The moving box is aligned
with Cartesian grid faces, which makes the integral evaluation straightforward
in an immersed boundary (IB) framework. Discontinuous and noisy derivatives of
velocity and pressure at the fluid-structure interface are avoided and
far-field (smooth) velocity and pressure information is used. We re-visit the
approach to compute hydrodynamic forces and torques through force/torque
balance equation in a Lagrangian frame that some of us took in a prior work
(Bhalla et al., J Comp Phys, 2013). We prove the equivalence of the two
approaches for IB methods, thanks to the use of Peskin's delta functions. Both
approaches are able to suppress spurious force oscillations and are in
excellent agreement, as expected theoretically. Test cases ranging from Stokes
to high Reynolds number regimes are considered. We discuss regridding issues
for the moving CV method in an adaptive mesh refinement (AMR) context. The
proposed moving CV method is not limited to a specific IB method and can also
be used, for example, with embedded boundary methods
Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library
In this paper we employ two implementations of the fictitious domain (FD)
method to simulate water-entry and water-exit problems and demonstrate their
ability to simulate practical marine engineering problems. In FD methods, the
fluid momentum equation is extended within the solid domain using an additional
body force that constrains the structure velocity to be that of a rigid body.
Using this formulation, a single set of equations is solved over the entire
computational domain. The constraint force is calculated in two distinct ways:
one using an Eulerian-Lagrangian framework of the immersed boundary (IB) method
and another using a fully-Eulerian approach of the Brinkman penalization (BP)
method. Both FSI strategies use the same multiphase flow algorithm that solves
the discrete incompressible Navier-Stokes system in conservative form. A
consistent transport scheme is employed to advect mass and momentum in the
domain, which ensures numerical stability of high density ratio multiphase
flows involved in practical marine engineering applications. Example cases of a
free falling wedge (straight and inclined) and cylinder are simulated, and the
numerical results are compared against benchmark cases in literature.Comment: The current paper builds on arXiv:1901.07892 and re-explains some
parts of it for the reader's convenienc
An effective preconditioning strategy for volume penalized incompressible/low Mach multiphase flow solvers
The volume penalization (VP) or the Brinkman penalization (BP) method is a
diffuse interface method for simulating multiphase fluid-structure interaction
(FSI) problems in ocean engineering and/or phase change problems in thermal
sciences. The method relies on a penalty factor (which is inversely related to
body's permeability ) that must be large to enforce rigid body velocity
in the solid domain. When the penalty factor is large, the discrete system of
equations becomes stiff and difficult to solve numerically. In this paper, we
propose a projection method-based preconditioning strategy for solving volume
penalized (VP) incompressible and low-Mach Navier-Stokes equations. The
projection preconditioner enables the monolithic solution of the coupled
velocity-pressure system in both single phase and multiphase flow settings. In
this approach, the penalty force is treated implicitly, which is allowed to
take arbitrary large values without affecting the solver's convergence rate or
causing numerical stiffness/instability. It is made possible by including the
penalty term in the pressure Poisson equation. Solver scalability under grid
refinement is demonstrated. A manufactured solution in a single phase setting
is used to determine the spatial accuracy of the penalized solution.
Second-order pointwise accuracy is achieved for both velocity and pressure
solutions. Two multiphase fluid-structure interaction (FSI) problems from the
ocean engineering literature are also simulated to evaluate the solver's
robustness and performance. The proposed solver allows us to investigate the
effect of on the motion of the contact line over the surface of the
immersed body. It also allows us to investigate the dynamics of the free
surface of a solidifying meta
An Immersed Interface Method for Discrete Surfaces
Fluid-structure systems occur in a range of scientific and engineering
applications. The immersed boundary(IB) method is a widely recognized and
effective modeling paradigm for simulating fluid-structure interaction(FSI) in
such systems, but a difficulty of the IB formulation is that the pressure and
viscous stress are generally discontinuous at the interface. The conventional
IB method regularizes these discontinuities, which typically yields low-order
accuracy at these interfaces. The immersed interface method(IIM) is an IB-like
approach to FSI that sharply imposes stress jump conditions, enabling
higher-order accuracy, but prior applications of the IIM have been largely
restricted to methods that rely on smooth representations of the interface
geometry. This paper introduces an IIM that uses only a C0 representation of
the interface,such as those provided by standard nodal Lagrangian FE methods.
Verification examples for models with prescribed motion demonstrate that the
method sharply resolves stress discontinuities along the IB while avoiding the
need for analytic information of the interface geometry. We demonstrate that
only the lowest-order jump conditions for the pressure and velocity gradient
are required to realize global 2nd-order accuracy. Specifically,we show
2nd-order global convergence rate along with nearly 2nd-order local convergence
in the Eulerian velocity, and between 1st-and 2nd-order global convergence
rates along with 1st-order local convergence for the Eulerian pressure. We also
show 2nd-order local convergence in the interfacial displacement and velocity
along with 1st-order local convergence in the fluid traction. As a
demonstration of the method's ability to tackle complex geometries,this
approach is also used to simulate flow in an anatomical model of the inferior
vena cava.Comment: - Added a non-axisymmetric example (flow within eccentric rotating
cylinder in Sec. 4.3) - Added a more in-depth analysis and comparison with a
body-fitted approach for the application in Sec. 4.
A fully resolved active musculo-mechanical model for esophageal transport
Esophageal transport is a physiological process that mechanically transports
an ingested food bolus from the pharynx to the stomach via the esophagus, a
multi-layered muscular tube. This process involves interactions between the
bolus, the esophagus, and the neurally coordinated activation of the esophageal
muscles. In this work, we use an immersed boundary (IB) approach to simulate
peristaltic transport in the esophagus. The bolus is treated as a viscous fluid
that is actively transported by the muscular esophagus, which is modeled as an
actively contracting, fiber-reinforced tube. A simplified version of our model
is verified by comparison to an analytic solution to the tube dilation problem.
Three different complex models of the multi-layered esophagus, which differ in
their activation patterns and the layouts of the mucosal layers, are then
extensively tested. To our knowledge, these simulations are the first of their
kind to incorporate the bolus, the multi-layered esophagus tube, and muscle
activation into an integrated model. Consistent with experimental observations,
our simulations capture the pressure peak generated by the muscle activation
pulse that travels along the bolus tail. These fully resolved simulations
provide new insights into roles of the mucosal layers during bolus transport.
In addition, the information on pressure and the kinematics of the esophageal
wall due to the coordination of muscle activation is provided, which may help
relate clinical data from manometry and ultrasound images to the underlying
esophageal motor function
- …